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3.2.1  Instability 


ABSTRACT 


The  propagation  of  rf  waves  launched  in  the  end  cell  and  central  cell  of  the  tandem  mirror 
propulsion  device  has  been  investigated  both  theoretically  and  experimentally.  Theoretically, 
a  computer  code  has  been  developed  to  study  the  wave  propagation  in  a  nonhomogeneous 
magnetic  Held.  It  was  found  that  the  amplitude  of  the  wave  excited  in  the  plasma  peaked 
while  approaching  the  resonance,  but  then  damped  out,  indicating  strong  absorption  of  the 
wave  by  the  plasma.  The  absorption  took  place  near  the  axis  and  midplane  of  the  device.  The 
experimental  results  confirmed  the  theoretical  prediction  of  the  phenomena  of  the  resonance 
effect  This  means  that  the  rf  power  is  heating  the  plasma  in  center  contrary  to  the  earlier 
prediction  that  the  heating  was  near  the  edge.  Therefore  higher  heating  efficiency  can  be 
possible.  A  very  imj>ortant  discovery  of  this  experiment  was  the  broadening  of  the  ICRF 
Fourier  spectrum  in  the  presence  of  the  plasma.  This  broadening  can  be  viewed  as  the 
detuning  of  the  Q  value  of  the  resonance  circuit  due  to  plasma  loading.  Instability  has  been 
observed  during  the  initial  phase  of  the  discharge  followed  by  quiescent  phase.  Therefore 
the  instability  was  stabilized.  Large  amount  of  data  has  been  taken  and  is  in  the  process  of 
being  analyzed.  The  analysis  of  plasma-gas  interaction  with  coaxial  injection  of  plasma  and 
neutral  gas  has  been  completed  which  will  he  reported  in  detail  separately.  A  collaboration 
effort  with  AVCO  Research  Laboratory  has  been  established  to  examine  the  feasibility  and 
method  of  using  the  high  temperature  plasma  of  this  device  for  reentrauce  study  for  national 
space  plane  developmental  effort. 


1  Introduction 

A  tandem  mirror  plasma  experimental  facility  was  established  for  space  propulsion  develop¬ 
mental  studies  1 ',  2;. 131,^1],  5j.  The  experimental  program  encompasses  three  areas  of  study. 
First,  to  establish  the  fundamental  plasma  properties  for  the  range  of  propulsion  application. 
Second,  to  seek  ways  to  maximize  plasma  heating,  thus  the  improvement  of  power  conversion 
efficiency  and  third,  to  study  the  plasma  propci lies  at  the  exhaust.  The  power  conversion 
method  will  be  to  inject  the  rf  power  Into  the  plasma  to  heat  the  ion  species  at  the  ion 
cyclotron  resonance  frequency  (ICRF).  Hence,  understanding  the  physics  of  wave  propaga¬ 
tion  in  this  device  is  fundamental  to  achieving  the  areas  just  mentioned.  Wa vq  propagation 
in  a  hot  plasma  contained  in  highly  nonhomogeneous  magnetic  field  is  a  very  complicated 
problem  which  was  hi  icily  touched  upon  during  the  closing  of  the  mirror  fusion  program. 
Nevertheless,  in  the  theoretical  area,  extensive  analytical  work  was  done  during  early  years 
of  this  program  leading  to  the  development  of  a  computer  code,  BEACH.  The  achievements 
in  theory  and  experiment  are  discussed  in  sections  2,  3,  *1  and  5.  The  industrial  participation 


and  the  proposed  scope  of  work  are  discussed  in  sections  5  and- 8. 


2  Progress  In  Theorectical  program 

2.1  Analytical  solution 

The  wave  propagation  in  a  cylindrical  geometry  for  cold  plasma  in  nonhomogencous  magnetic 
field  has  been  solved  analytically  and  numerically.  The  mathematics  involved  Is  long  and 
complex,  berefore  the  details  will  be  discussed  in  a  separate  report  rcfjGj.  Here,  the  greatly 
reduced  but  yet  comprehensible  theory  is  presented. 

It  can  be  shown  that  in  the  vacuum  the  Maxwell's  equations  can  he  reduced  to  the 
following  differential  equation: 


The  solutions  to  the  above  differential  equations  are  the  Messel  functions  ln[vr)  and  h'n(vr). 

It  is  seen  that  Maxwell's  equations  can  he  separated  into  two  component  solutions  in 
a  vacuum:  TIv  and  TM.  While  the  presence  of  a  plasma  prevents  a  complete  separation 
into  the  two  modes,  the  Maxwell’s  equations  can  he  simplified  into  two  similar  modes:  the 
“nearly”  TIv  mode, 


and  the  electric  nmdesfnearly  TM), 
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The  plasma  field  solutions  are  then  simply  a  sum  of  the  nearly  TE  and  nearly  TM  modes. 
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Where  the  solution  of  Equations  3  and  6  is  the  Bessel  function  ./„ 

lUx  -  HrlnM 
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Given  the  boundary  conditions  and  the  transverse  fields,  the  solution  can  be  found  to  the 
power  absorption  problem.  However,  the  form  of  the  solution  is  extremely  complex  and  any 
physical  insight  is  lost  among  the  algebra.  In  order  to  provide  some  physical  insight  into 
the  solution,  it  is  necessary  to  reduce  the  problem  further.  Thus,  several  more  assumptions 
are  mftde‘0)E*  ~  Q,i.e.,onlv  TE.  modes  are  considered;  (2)  n  -  0,  axisymmetric  modes;  (3) 
nK(r)  iij(r)  Constant.  Therefore  the  fields  in  the  plasma  are  reduced  to  the  following 
form: 
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where  the  variables  A.l>.  and  S  are  defined  in  ref  li .  The  absorbed  plasma  power  was  shown 
to  be 
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2.2  Computational  Analysis 

Although  the  analytical  treatment  of  the  previous  section  yielded  intuitive  results  of  the 
dependence  of  Paf,s  on  plasma  parameters,  it  nevertheless  is  loo  simplified  a  treatment  to 
apply  to  a  real  system.  The  main  simplification  that  was  made  in  the  analytical  treatment 
was  the  use  of  only  one  eigenmode.  Although  that  allowed  for  a  simple  solution,  the  single 
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mode  treatment  is  nevertheless  too  crude  to  be  used.  A  typical  solution  that  proved  to 
be  useful  required  the  solution  of  about  30  50  modes  which,  upon  summation,  form  the 
resulting  field  solutions. 

Another  limitation  of  the  analytical  appinach  is  that  it  doesn’t  consider  the  issue  of  the 
field  inhomogeneity.  Again,  this  simplification  was  made  in  order  to  reduce  the  complexity  of 
the  problem.  It  is  thus  for  these  reasons  that  a  computational  analysis  was  deemed  necessary 
and  undertaken. 


The  goal  of  this  computational  analysis  is  to  provide  a  better  modeling  of  the  wave 
plasma  coupling  problem  than  the  simplified  results  of  the  analytical  solution.  However,  this 
computational  analysis  is  not  simply  a  purely  numerical  simulation.  Much  analytical  work 
is  done  to  reduce  the  fields  to  functional  form  ,  i.e.  ln[kTr)%  h’tx[krr ),  J,x(ki\r)  before  the 
numerical  "number-crunching”  is  handed  over  to  the  codes. 

The  BEACH  code  extends  the  analytical  approach  taken  in  the  cylindrical  geometry  so¬ 
lution  of  the  previous  section.  It  solves  for  the  IC'RP  wave  propagation  in  an  inhomogeneous 
axial  magnetic  field.  The  inhomogeneous  axial  geometry  is  handled  by  separating  the  axial 
length  into  discrete  axial  “slices",  each  of  which  contains  a  homogeneous  axial  magnetic  f  eld. 

The  code  first  finds  the  eigenmodcs  that  exist  in  each  axial  slice  by  satisfying  both  the 
infinite  plasma  dispersion  relation  and  the  boundary  conditions.  The  boundary  condition  of 
the  plasma  vacuum  system  can  be  written  as: 


The  solution  to  the  above  are: 


C\IIp  t  Cj  /#j, 

(25) 

(20) 

"v 

;<?!<•,  r.'ivj  , 

/»  '>■* 

(27) 

// 

(28) 

D 
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(22) 

The  term  D  is  then  the  boundary  condition  of  the  bounded  plasma  vacuum  system.  Solving 
this  boundary  condition  coupled  with  the  dispersion  relation 

Sk[  +  k]  | k](S  +  P)  ~  k*(SP  4  RL)\  +  P[k\  -  klR){kl  -  klL)  -  0  (30) 
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yields  the  eigen  inodes  of  the  cylindrical  system. 

After  solving  for  the  eigenmodcs  of  the  system,  the  code  then  computes  the  excitation 
coefficients  of  the  initial  axial  slice  as  driven  by  the  source  term  (Antenna).  The  excitation 
coefficient  of  the  next  axial  slice  can  'hen  he  determined  by  imposing  two  conditions  of 
continuity  on  the  boundary  between  the  two  axial  slices:  conservation  of  the  integrated 


axially  propagating  Poynting  Flux  f  !>tr<lr  ,  and  the  continuity  of  wave  phase.  These  two 


conditions  uniquely  determine  the  excitation  coefficients  of  the  different  eigenmodcs  in  the 


next  axial  slice.  This  analysis  is  then  performed  for  all  axial  slices. 


2.3  Computational  Results 

Computational  analysis  was  performed  on  two  types  of  plasma  regimes:  collisionless  and 
cullisional  cold  plasma.  Different  physics  arc  inherent  in  each  different  regime. 


2.3.1  Collisionlcss 


Initially,  the  calculations  of  the  ICTlF  field  propagation  were  performed  using  a  cold,  colli- 
sionless  plasma  model.  This  ensures  that  there  exist  no  complex  k .  roots  that  could  damp 
out  the  axially  propagating  Poynting  Flux,  thus  enabling  a  good  check  of  the  model's  ability 
to  conserve  power. 


Several  different  antenna  geometries  are  used  to  launch  different  modes.  The  two  antenna 
geometries  that  were  studied  were  the  full  turn  loop  and  the  half  turn  loop. 

I'he  baseline  mode)  is  driven  by  a  Full  Turn  lamp  ( FTI.)  antenna.  Because  it  is  symmetric 
in  the  full  turn  loop  antenna  excites  only  the  n  0  mode.  Figures  1  and  2  show  the 
propagation  characteristic  *»f  the  K'llF  IJ.  and  E}  fields  as  the  waves  propagate  from  the 
launch  point  (/.  (J)  towards  the  resonance  point  downstream.  The  results  that  are  readily 
observable  from  ibis  collisionless  propagation  of  the  n  U  mode  towards  resonance  arc  given 
bel«  rtv: 


1 .  Ii.  Increases  -  This  is  an  expected  result.  Stix  shows  (Stix,  Theory  of  Plasma  Waves, 
p.51.  Eq(18)  j  H  that  for  a  homogeneous  lossless  r.icdium,  the  group  velocity  is 


N=0,  N&=  1.50x  JO1',  >/WClO=  0.70$  ,  COLD  COLUSIONLKSS  MODEL 
NO.  MODES  =  30 


Figure  1:  IC’RF  B*  Propagation  in  a  Cold  Collisionless  Plasma,  Full  Turn  Loop  Antenna, 
u.'  u.'eiU  ■■  U.788.  Waves  are  launched  by  the  ICRF  antcnm  at  r  =  Ocm  and  propagate 
towards  resonance  at  z  -  50cm. 
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N=0.  NEX=  1.50hJ0,,(  TV/WCIO^  0700  ,  COLD  COLLISIONLESS  MODEL 
NO,  MODES  -  *30 


Figure  2:  IC'liF  E %  Propagation  in  a  Cold  Collisionless  Plasma,  Full  Turn  Loup  Antenna, 
-  0.788 
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When  resonance  is  approached,  the  group  velocity  slows  down,  resulting  in  an  increased 
cnergv  density  U*„,  This  increased  energy  density  is  apparent  in  the  increase  in  1/i.j, 
since  the  energy  density  varies  as 


...  /*:*  /ir 

"  w  /'o  *  <0  () 


2.  Increase  in  axial  and  radial  k*  ■  As  resonance  (w  u/„)  is  approached,  the  left 
hand  cut-off  (A*?  k^L)  increases,  causing  an  increase  in  the  resulting  A\  cigcnmodcs. 
This  is  consistent  with  the  slowing  down  of  the  group  and  phase  velocities  as  resonance 
is  approached. 


3.  High  edge  electric  fields  Et ,  E  and  E »  that  propagate  inward  •  This  is  a 
result  that  could  be  consistent  with  experimental  observation.  Golovato  points  otu 
(Persona)  communications)  that  although  a  high  edge  electric  field  is  predicted  by  the 
ANTENNA  code,  experimentally  no  edge  electric  field  is  detectable.  The  BEACH 
simulation  shows  that  the  reason  that  no  edge  electric  field  is  detectable  is  that  the 
edge  spike  moves  radially  inward  into  the  plasma  as  resonance  is  approached. 

•1.  Higher  order  modes  excited  at  higher  amplitudes  for  lower  initial  u>/u/ri  -  This 
result  reflects  the  effect  of  the  antenna  excitation  profile  ( Jt(ks ))  on  the  excited  EM 
field  coefficients.  For  a  lower  w/a,*,.,,  the  value  of  the  left-handed  cut*ofr  in  the  infinite 
plasma  dispersion  relation  A*/,  decreases.  Since  the  ./#(A\)  profile  remains  constant, 
this  means  that  ALL  modes  are  now  excited  to  higher  values.  (See  Figure  ??  ) 


I  he  solution  of  the  dual  half  turn  loop  (DHTL)  geometry  puses  a  more  difficult  problem 
than  the  full  turn  loop  (F  PL)  solution.  'Die  main  difference  between  the  UHTL  and*  the 
FTL  lies  in  the  ./(A*.)  driving  term,  specifically  the  radial  feeders.  This  is  different  from  the 
Full  Turn  Loop  (FTL)  geometry  in  which  the  contribution  from  the  radial  feeders  cancel 


out  (Since  one  leg  is  entering  and  the  other  is  exiting,  net  Jr  is  zero).  The  modelling  of  the 
driving  term  J{k;)  for  the  DIITL  is  detailed  in  a  full  report  in  preparation. 

The  basic  characteristics  of  ICRF  propagation  towards  resonance  are  similar  to  the  FTL 
(n-0)  case,  i.e.  increasing  A*,  and  kx.  and  increasing  3 .  amplitude.  The  main  difference 
is  the  larger  Av  and  A*.  ior  the  DHTL  case,  which  is  to  be  expected  since  the  asymmetric 
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(about  t  0)  DIITL  excitation  necessitates  radial  eigenmodcs  that  have  smaller  perpendicular 
wavelengths  (larger  A*a)  and  hence  larger  A\. 

The  £i  Held  also  appears  to  he  similar  in  structure  to  the  FTL  case,  with  an  initial  edge 
peaked  that  propagates  towards  the  center  of  the  plasma  as  resonance  is  approached. 


2.3.2  Collisioual 


The  presence  of  collisions,  of  course,  implies  a  Unite  temperature.  The  Cold  Pl»sim  Colli- 
sionnl  Model  means  that  the  expansion  of  the  plasma  Z  function  excludes  temperature  ejects 
and  keeps  only  terms  of  t)(i/,  *;).  The  code  results  for  the  case  of  ICRF  propagation  in  a 
cold  collisioual  plasma  is  seen  below  in  Figure  3. 

The  presence  of  collisions  manifests  itself  in  a  complex  dielectric  tensor.  Recall  that  the 


dielectric  tensor  is  given  by 


where 


S  W  0 

i!)  S  0 

0  0  r 


] )  wr'w"  (*J  21 J 

".T  ^ 


V  WJ" 

—  tl  rll'.ol 


(37) 

(38) 


Using  these  complex  elements  of  the  dielectric  lensor,  the  dispersion  relation  can  be  solved 
with  the  geometric  boundary  condition  to  yield  a  set  of  cigenmode  solutions  for  this  plasma 
vacuum  geometry.  The  difference  between  these  solutions  and  those  of  the  cold  collisionless 
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N=0,  N£=  1.00x30'*,  W/WCiO=  0.780  ,  COLD  COLL1SIONAL  MODEL 
T)-PERP=  10.00,  TE-PERP=  10.00,  Tl~Z=  10.00,  TE-Z-  10.00 
NO.  MOPES  =  30 


plasma  is  that  the  k.  cigeumodcs  are  now  comph  x  these  waves  jure  now  damjKtd  as  they 
propagate  towards  the  cyclotron  resonance.  This  damping  transfers  energy  from  the  waves 
to  the  plasma,  causing  plasma  heating. 


3  Progress  In  Experimental  Program 

3.1  Diagnostics 

hr  order  to  determine  the  plasma  characteristics  of  the  discharge  produced,  an  adequate  set  of 
diagnostics  must  he  assembled  to  measure  some  basic  plasma  parameters.  These  parameters 
include  th  a  temperature  i‘„  electron  temperature  Tc,  and  electron  and  ion  densities  n,, 
nr.  A  system  of  diagnostics  was  designed,  fabricated  or  reconstructed,  calibrated  and  used 
to  determine  plasma  properties  as  part  of  this  experiment.  This  section  will  disemut  the 
diagnostics  systems  that  were  fabricated  and/or  reconstructed  for  use  on  this  experiment. 

The  diagnostics  currently  operating  are  listed  in  Table  2  along  with  the  plasma  properties 
that  they  measure.  The  principle  of  these  diagnostics  and  their  operation  and  mrtfcod  of 
analysis  are  evaluated  and  presented  below. 


Table  2:  jPPEX  Diagnostic  Set 


DIAGNOSTIC’ 


.  LOCATION  |  PLASMA  PROPERTY  MEASURED 


LANGMUIR  PROBE  CC  Tr(u J.Mu) 

IN  I  KRh'KRo.MK  I  KR  CC.  Anchor  nf 
//„  ARRAY  CC.  Anchor  n(r) 

SPECTROMETER  ,  Anchor  j  7',,  Plasma  Rotation 

DIAMAGNETIC  LOOP  ;  CC  -  J.  -  Te 

i 

j  13-DOT  PROBE  j  C’C 

1  i 


k  ,  IC’RF  Propagation 


3.1.1  Langmuir  Probe 


The  basic  principle  of  the  Langmuir  is  well-known  and  is  briefly  described  to  the  method 
of  analysis  used  here.  A  Current- Voltage  characteristic  curve  k  produced  when  the  probe 
driving  voltage  is  sweeped  from  some  negative  value  to  a  positive  value  (ilk  b,  p.ll5)refj9j. 
At  negative  voltages,  the  probe  should  attract  only  ions  at  an  ion  saturation  flux  found  to 
be  (Huddleston.  Leonard,  p.150)  iteffl 
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where  A  is  the  probe  area,  «,<  the  plasma  density,  and  M  the  ion  mass.  As  the  probe  voltage 
is  increased,  more  and  more  electrons  can  penetrate  the  negative  potential  until  a  floating 
potential  \)  is  eventually  reached.  The  probe  being  at  the  floating  potential  means  that 
the  net  inward  electron  flux  equals  the  net  inward  ion  flux.  As  the  voltage  is  continuously 
raised,  the  probe  now  attains  a  linear  regime  in  which  the  electron  current  behaves  in  an 
exponential  manner,  assuming  a  Maxwellian  electron  distribution  (Ilk  L,  p.J35).  This  is 
known  as  the  transition  region  and  has  the  characteristic: 


/ 


.  (  kT  \  -  .* 

■V'“  (,*,„)  <"• 


From  experimental  observations  of  the  current -volt age  characteristic,  the  electron  tem¬ 
perature  l'r  and  density  ne  can  be  found.  To  calculate  the  electron  temperature,  the  slope  of 
the  transition  region  is  used.  ‘Inking  the  region  of  constant  slope  in  this  region,  the  electron 
temperature  can  be  found  from  (Hutchinson,  p.0-1)  ref  HI 
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where  1  is  the  Langmuir  probe  current  in  the  middle  of  the  linear  region,  is  the  slope  td 
tile  Langmuir  Probe  1  V  characteristic  of  the  linear  region,  and  ls,  the  ion  saturation  current. 

From  'ids  value  of  the  electron  temperature,  the  edge  plasma  density  can  be  found  from 
the  ion  saturation  current  Jit  (Hutchinson.  p.Utl): 
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where  .1,  is  the  area  of  the  Langmuir  Probe  tip,  C\  the  ion  sound  speed  (0,  (T,/w,)’  ). 


Tire  densitv  is  then  found  to  be 


••UV/ 


3.1.2  Interferometer 


The  basic  principle  of  plasma  interferometry  utilizes  the  phase  shifting  of  the  plasma  CLmodc 
wave  whose  <lispersion  relation  is 


J  JL  *  W 


where 


is  the  electron  plasma  frequency.  The  index  of  refraction  can  then  be  defined  as  : 

cU  up  * 

A*  -  1  -  " 

The  phase  change  due  to  the  wave  propagating  in  the  plasma  medium  of  length  l  (in  addition 
to  its  normal  propagation  in  vacuum)  is  then 

Art  kl(  .V  1) 

Prom  the  measured  value  of  the  A$,  the  line  averaged  electron  density  can  be  found  by; 
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3.1.3  Frequency  Sweeping  Operation 


The  above  is  fur  a  homodym*  interferometer.  The  interferometer  used  in  the  PPEX  experi¬ 
ment  is  a  frequency-swept  interferometer.  The  basic  principle  of  the  phase  shift  is  the  same 
as  that  of  the  homodvne  interferometer,  the  difference  being  that  the  frequency  of  the  in- 
terferometer  is  swept  over  a  fixed  range.  Sweeping  frequencies  from  HIUKHz  to  l.MIJz  have 
been  used.  The  interferometer  in  the  PPEX  experiment  utilizes  a  SOKliz  ramp  generator  to 
produce  the  sweeping  frequency. 


1G 


For  such  a  frequency  swept  interferometer,  the  phase  shift  is  given  by 
A#  2a-  x  1.3*1  x  10  ,tf7yi,r  .  /«•(*» i)rf*(w'a) 

For  the  35GHz  interferometer  used  on  the  PPEX,  this  works  out  to  be 

Atf  3..X1  x  10  *V Ax 
or,  using  an  average  plasma  radius  of  10cm, 

/»,  A#  «  (2.UH  x  ItFViu  3)  (13) 

where  A#  is  the  phase  shift  in  radians. 


3.1.4  Spectrometer 

The  Doppler  shift  due  to  a  source  moving  parallel  (with  velocity  vj  to  the  direction  of 


observation  is 


AA  : V’  A 


It  can  be  shown  (H&1..  p.2f>8)reP9  that  for  sources  of  Maxwellian  velocity  distribution, 


/pA)  AM  ut,). 


where  It  is  the  total  line  intensity.  The  intensity  is  half  its  maximum  value  when  the  expo* 
nential  reaches  ’.  The  full  width  half  max  (l*\Vli.\t)  (llJtr  1..  p.'JGtl)  line  width  is  found  to 


A.\i  /(  T.1U  10  7\[T'M)j 

(/,  is  in  degrees  Kelvin,  M  is  the  molecular  wt  of  the  ion)  Therefore  the  ion  temperature  is 


found  to  be 


AA^^ 
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The  rotating  mirror  assembly  is  constructed  to  convert  the  J*A  Monochromcter  Spectrom¬ 
eter  into  a  scanning  spectrometer.  The  reason  for  using  a  rotating  mirror  assembly  is  to 
translate  the  incoming  beam  by  a  glass  slide,  hence  shifting  the  -/avclcngth  entering  the 
photomultiplier  slit. 


Tsing  Snell’s  Law,  the  amount  of  beam  translation  is  found  to  be 


um  h  ,  ™e'-  ) 

\  \'1'  (sinfi.JV 

By  translating  the  beam  with  a  rotating  glass  slide,  the  wavelength  that  is  being  examined 
by  the  photomultiplier  assembly  is  thus  caused  to  vary. 


3.1.5  //„  Army 

The  basic  //„  array  constructed  for  use  on  the  IM’KX  is  a  multi-chord  photomultiplier  system. 
The  basic  principle  of  the  Muliichord  Array  is  that  of  a  pinhole  camera.  Light  from  the 
plasma  passes  through  a  pinhole  and  is  subsequently  detected  by  the  photodetector  array 
on  the  end  of  the  array  holder.  Arbitrary  filters  may  be  installed  over  the  pinhole  to  allow 
viewing  of  different  light,  he  it  impurity  or  plasma  emissions. 

3.1.0  Diamagnetic  Loop 

The  diamagnetic  loop  is  a  magnetic  field  probe  that  is  positioned  around  the  plasma.  The 
loop  detects  the  magnetic  flux  that  it  surrounds  and  the  voltage  induced  on  the  loop  is 
subsequently  measured  and  processed  to  yield  the  plasma  beta. 

Tsing  Maxwell’s  equations  and  the  MHO  Equilibrium  condition,  it  can  be  shown  that 
the  magnetic  field  inside  a  plasma  of  finite  J  is  given  by 

/i,  v'1  Jli» 

where  J  is  defined  as  the  ratio  of  the  plasma  pressure  to  the  magnetic  field  pressure,  J 
I'  (If  The  llnx  that  passes  through  the  diamagnetic  loop  is  then  found  to  be 


O.  /<•(«(  .V 


The  plasma  eontribulion  is  then  seen  to  be 


.XJU*-r* 


The  output  voltage  that  is  induced  in  the  loop  is  RC  integrated  to  yield  the  magnetic  flux 


and  is  seen  to  lie 


i  >  $*.1’ 

,0  lie 
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Thus  Ur*  plasma  0  can  he  found  directly  from  this  output  voltage  and  is  seen  to  be 

’IRC 

,,  i  • 

V  »}2^  jTl! 

A  B&wr l 

where  the  factor  -f  takes  into  account  the  frequency  response  of  the  diamagnetic  h>op. 


3.1.7  B-Dol  Probe 

The  B»Dot  probe  picks  up  the  IC’HF  fluctuating  magnetic  Helds  that  are  either  propagating 
in  the  plasma  or  evanescent  but  still  delectable.  The  13-Dot  I’robe  works  with  a  similar 
principle  to  the  Diamagnetic  Loop,  namely  that  the  lime-varying  RF  field  is  picked  up  and 
read  as  an  induced  loop  voltage  V,  where 

This  induced  loop  voltage  can  then  be  detected  and  digitized. 

3.2  ICRH  Wave  Experimental 

3.2.1  experiment  set-up 

The  schematic  of  the  experiment  is  shown  in  Fig.  5.  There  are  double  half-loop  antennae  in 
the  south  end  cell  and  in  the  north  end  of  the  central  cell.  Doth  antennae  are  located  at  fields 
higher  than  the  corresponding  ones  f«r  their  resonant  frequencies  (a,*  u/n).  The  frequency 

uf  the  rf  power  is  in  resonance  with  the  iun  cyclotron  frequency  at  the  field  of  the  midplane 
of  each  cell.  When  rf  radiation  was  emitted  at  the  antennae,  the  wave  excited  in  the  plasma 
travels  both  radially  inward  and  axially  toward  the  center  of  each  cell.  The  diagnostics 
enclosed  by  the  box  in  Fig.  o  are  currently  operational  which  were  built  and  installed  in 
the  past  year.  A  scanning  spectrometer  and  a  retarding  field  grid  energy  analyzer  will  be 
installed  later. 


3.2.2  experimental  results 

The  B-dot  probe  allows  the  measurement  of  /)r,  IJg  and  B..  13r  and  Bg  were  very  small. 
The  Fourier  spectrum  of  B .  is  shown  in  Fig.  0  is  the  wave  in  the  vacuum.  The  fundamental 
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frequency  is  3  MHx  marked  by  II.  The  harmonies  arc  marker)  by  1.2,3  am)  4.  The  jnraks 
in  between  may  be  due  to  the  imperfection  of  the  rf  transmitter.  A  comparison  of  Figs.  6a 
and  6b  shows  that  the  amplitude  of  fundamental  frequency  peak  U  is  one  unit  in  the  vacuum 
and  50  units  in  the  plasma.  The  half-width  is  about  0.1  MHx  for  the  vacuum  and  1  Mil* 
for  the  plasma.  Therefore  the  amplitude  of  the  wave  is  50  times  higher  in  the  phnma  than 
in  the  vacuum  at  the  resonance  which  confirms  the  theoretical  prediction.  The  width  of 
the  wave  is  10  times  broader.  The  reason  that  the  plasma  signals  are  much  larger  than 
the  vacuum  signals  is  that  the  ICRF  IJA  is  evanescent  in  the  vacuum,  while  the  presence 
of  the  plasma  allows  it  to  propagate.  This  large  broadening  indicates  increasing  resistance 
due  to  the  loading  of  plasma.  There  is  little  evidence  that  the  wave  injected  in  the  end  cell 
will  propagate  into  the  central  cell  or  vice  versa  because  of  the  fact  that  the  wave  energy  is 
damped  and  absorbed  near  the  resonance  in  the  respective  cells.  There  is  also  little  evidence 
of  heating  electrons  at  the  edge,  an  indication  of  less  loss  ami  higher  efficiency. 


3.2.3  Plasma  Broadening  of  ICRF  Spectrum 


One  of  the  discoveries  made  on  this  plasma  propulsion  experiment  is  the  broadening  of  the 
ICRF  Fourier  spectrum  in  the  presence  of  a  plasma.  Comparison  of  a  vacuum  ICRF  shot 
to  a  plasma  ICRF  shot  clearly  shows  this  effect.  We  propose  an  explanation  for  this  effect. 
The  broadening  is  thought  to  he  caused  by  the  detuning  of  the  plasma  Q  sis  defined  by 


/  -• 
A/  R 


Sime  the  atiiemm  was  originally  tuned  for  vacuum,  the  Q  is  very  large  for  si  vsicuum  shot 
(hence  a  narrow  peak).  However,  the  presence  of  a  plasma  increases  the  serial  resistance 
of  the  equivalent  RFC  resonance  circuit,  therefore  Q  decreases,  causing  A/(  to 
increase. 


3.2.4  Instability 

To  understand  the  resonance  effect  two  types  of  experiments  were  carried  out:  To  inject  the 
rf  power  sit  a  frequency  slightly  higher  or  lower  than  the  resonance  frequency.  Therefore 
we  'ailed  the  class  of  discharges  with  u>  >  u-v,  as  type  1  and  the  class  of  discharges  with 
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Fourier  Power  Amplitude 


FREQUENCY  (MHz) 


FREQUENCY  (MHz) 

Figure  9:  Fourier  spectrum  of  the  measured  wave  of  frequecy  of  3  A/ if,:  (a)  m  the  vacuum  and 

(b)  in  the  plasma.  The  Amplitude  and  width  of  the  wave  in  the  plttm*  are  mrh 
Ivfer  then  in  the  vacuum. 


u?  %  u?n  its  type  II.  Figure  7  shows  the  H-a  acquired  during  TYPE  I  and  11  discharges. 
There  are  distinguish  features  between  these  two  types  of  H-a  emissions.  For  the  first  type, 
There  is  strong  fluctuations  initially  which  decays  and  becomes  quiescent.  This  indicates  the 
instability  at  the  initial  phase  of  the  discharge.  Fortunately  it  was  stabilized  possibly  due 
to  the  ponderomotive  force  of  the  rf  wave.  The  plasma  is  quiescent  for  the  type  II  discharge 
although  the  amplitude  of  the  11-a  emission  is  modulated.  This  modulation  may  he  due  to 
the  ECRU  pulse  or  plasma  rotation  which  has  to  he  carefully  investigated.. 

T«  find  out  the  physical  nature  of  the  instability  and  its  stabilizing  mechanism  for  the 
type  l  discharge  and  source  of  the  modulation  for  the  type  11  discharge  will  be  the  center  of 
<>ur  research  effort.  T«  do  so  additional  diagnostics  are  needed,  such  as  13-dot  probe  array, 
Langmuir  probe  array,  retarding  field  grid  energy  analyzer  and  charge  exchange  analyzer. 

In  addition  to  the  results  presented  above  large  amount  of  data  have  been  collected. 
Detailed  analysis  is  underway. 
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H-oipho  Choroclerislic  for  Shot  No  1854 


590  595  bl<J  605  610  615  620  625  630  635  640 


Time  (msec) 


Figure  7:  Comparison  between  the  11-a  output  of  TYPE  1  and  TYPE  II  discharges.  TYPE  1 
(Top)  shows  a  much  higher  light  emission  level  as  well  as  the  three  different  regions  identified 
as  Ramp-up,  Oscillatory,  and  Quiescence.  TYPE  II  (Bottom)  shows  little  light  emission  and 
is  relatively  quieiccnt 
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